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I. Introduction 



The similarity assumption reduces the complexity of partial differential equa- 
tions as it turns the governing partial differential equations to relatively sim- 
ple ordinary differential equations. The Einstein field equations (EFEs) 

Rab - ^QabR ^ t^Tab, (o, 6 = 0, 1, 2, 3) (l) 

are non-linear partial differential equations. Self-similar solutions have been 
shown to be very useful in solving these set of equations. 

Self-similarity refers to the fact that the spatial distribution of character- 
istics of motion remains similar to itself at all times during the motion. Sim- 
ilarity solutions were first studied in General Relativity (GR) by Cahill and 
Taub [1] . They assumed that the solution was such that the dependent vari- 
ables were essentially functions of a single independent variable constructed 
as a dimensionless combination of the independent variables. In the simplest 
situation, a similarity solution is invariant under the transformation r ar, 
t ^ at for any constant a. Geometrically, they showed that the existence of 
a similarity of first kind in this situation could be invariantly formulated in 
terms of the existence of a homothetic vector. A natural generalization of 
homothety, called kinematic self-similarity, exists and is defined by the exis- 
tence of a kinematic self-similar (KSS) vector field. A KSS vector ^ satisfies 
the following conditions: 

£^hab = 25/la6, (2) 
£^Ua = aua, (3) 

where hat is the projection tensor, and a and 6 are constants. The similarity 
transformation is characterized by the scale- independent ratio a/ 6, which is 
known as the similarity index. 

By using a similarity index. Carter and Henriksen [2,3] defined other kinds 
of self-similarity: namely, second, zeroth, and infinite kinds. In the context 
of kinematic self-similarity, homothety is considered as the first kind. Sev- 
eral authors have explored KSS perfect-fiuid solutions. The only barotropic 
equation of state is compatible with self-similarity of the first kind is p = kp. 

Carr [4] has classified self-similar perfect-fiuid solutions of the first kind 
for the dust case {k = 0). The case < A; < 1 has been studied by Carr and 
Coley [5] . Coley [6] has shown that the Friedmann- Robertson Walker (FRW) 
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solution is the only spherically symmetric homothetic perfect-fluid solution 
in the parallel case. Mcintosh [7] has discussed a stiff ffuid {k = 1) being the 
only compatible perfect-ffuid with homothety in the orthogonal case. Benoit 
and Coley [8] have studied analytic spherically symmetric solutions of the 
EFE's coupled with a perfect-fluid and admitting a KSS vector of the flrst, 
second, or zeroth kind. 

Carr et al. [9] considered the KSS vector associated with the critical be- 
havior observed in the gravitational collapse of spherically symmetric perfect 
fluid with equation of state p = kp. Carr et al. [10] further investigated 
the solution space of self-similar spherically symmetric perfect-fluid models 
and the physical aspects of these solutions. They combined the state space 
description of the homothetic approach with the use of the physically in- 
teresting quantities arising in the co- moving approach. Maeda et al. [11] 
discussed the classiflcation of the spherically symmetric KSS perfect-fluid 
and dust solutions. Recently, Sharif and Sehar investigated the classiflcation 
of cylindrically symmetric [12] and plane symmetric [13] KSS perfect-fluid 
and dust solutions. 

The existence of self-similar solutions of the first kind is related to conser- 
vation laws and to the invariance of the problem with respect to the group of 
similarity transformations of quantities with independent dimensions. This 
can be characterized in GR by the existence of a homothetic vector. Per- 
veen [14] classified plane symmetric Lorentzian manifolds according to their 
homotheties and found different solutions admitting 5, 7, or 11 homotheties. 
Among these solutions, two correspond to 5 homotheties and five admit 7 
homotheties. The only solution admitting 11 homotheties is the Minkowski 
metric. 

Recently, Sharif and Sehar explored the physical properties of spherically 
symmetric self-similar solution of the first kind [15] and cylindrically sym- 
metric self-similar solution of the first kind [16]. This is the third paper 
in the series. Here, we are extending the same analysis for the plane sym- 
metric self-similar solution of the first kind. The paper can be outlined as 
follows: In Section II, we shall write down the self-similar solutions of the 
plane symmetric spacetime. Section 111 is devoted to a discussion of the 
physical properties of these solutions both in co-moving and non-co-moving 
coordinates. In section IV, we shall explore the singularity structure of these 
solutions. Finally, we shall summarize and discuss all the results in Section 
V. 
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II. Plane Symmetric Self-Similar Solutions of 
the First Kind 

The general plane symmetric spacetime is given by the line element [17] 

ds^ = e^-(t'-^dt^ - e^^^^'^^dx^ - e^''^''^\dy^ + dz''), (4) 

where u, A and /i are arbitrary functions of t and x. Perveen [14] classified 
plane symmetric Lorentzian manifolds by homotheties and found self-similar 
solutions of the first kind. There are two classes of such solutions, one ad- 
mitting 5 homotheties and the other admitting 7 homotheties. This paper is 
devoted to discussing the physical properties of these solutions. 
The first metric is given by 

ds^ = e'^^^^Uf - dx^ - e^'^^^^dy'' + dz^). (5) 

This metric has the following two solutions having 5 and 7 homotheties, 
respectively: 

ds^ = {^f^dt^ - dx^ - (-)'^(dy' + dz^), (6) 

Xo Xo 

where B, 5^0, ^ 7^ 1, and 

ds'' = {- f^{dt^ - dy^ - dz") - dx\ (7) 

Xo 

where Ay^O. The second metric has the form 

ds^ = df - e^-^Wdx^ - e2M*)(rf^2 ^ ^^2y 

This metric has the following two solutions with 5 and 7 homotheties, re- 
spectively: 

ds'' = dt^ - {-f^dx" - i-y^'idy^ + dz^), (9) 
where Ay^B, B 0, A ^ 1, and 

ds^^dt^-{-f''{dx'' + dy'' + dz''), (10) 

^0 

where A ^ 0. 
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The metric 

where a 7^ 0, has only one solution admitting 7 homotheties, and that solution 
is given by 

ds^ = i-fdt' - dx^ - {-fe^-a{dy' + dz^). (12) 

Xq Xq 

The metric 

ds^ = _ g2/W^^2 _ ^2^+2/(0(^^2 ^ ^^2)^ (13) 

where a ^ 0, also has one solution with 7 homotheties, and that solution is 
given by 

ds' = df - i-fdx' - i-fe^Hdy^ + dz"). (14) 
Finally, the metric 

ds" = de - dx^ - e2(»*+''^) (ciy2 + dz"), (15) 

yields the following solution with 7 homotheties: 

ds" = de - dx^ - e'<'+''\dy'' + dz"). (16) 

Thus, we have total of seven self-similar solutions that are given by 
Eqs.(6), (7), (9), (10), (12), (14), and (16). These can be divided into two 
classes admitting 5 and 7 homotheties, respectively 



III. Kinematics of the Velocity Field 

This section is devoted to discussing the kinematical properties of the self- 
similar solutions of the first kind both in co-moving and non-co-moving coor- 
dinates. These properties [17] can be listed as follows: The volume behavior 
of the fluid can be determined by the expansion scalar defined by 

e = u%. (17) 

The acceleration can be defined as 

Ua = Ma;fo«^ (18) 
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where Ua is the four-vector velocity. The rotation is given by 

l^ab = Ula-b] + U[aUb]. (19) 

The shear tensor, which provides the distortion arising in the fluid flow leav- 
ing the volume invariant, can be found by 

CTab = ^i{a;6) + U(^aUb) " ^©/iab- (20) 

The shear scalar, which gives the measure of anisotropy, is defined by 

a = (7„^.^7"^ (21) 

The expansion rate with respect to proper time is given by Raychaudhuri's 
equation [18] 

^ = -^e^ - aaba'''' + (^abU^u' - RabU^uK (22) 
otr o 

We now discuss these properties for the self-similar solutions given in the 
previous section. 



1. Kinematic Properties in Co- moving Coordinates 

First, we evaluate the kinematical properties of the self-similar solutions in 

co-moving coordinates. For this purpose, we have divided the solutions into 
two classes. Class 1 has solutions that admit 5 homotheties while class 2 has 
solutions that admit 7 homotheties. 



A. Class 1 

This class has two solutions given by Eqs.(6) and (9). We start the discussion 
of kinematic properties for the solution given by Eq.(6). 

The expansion scalar is zero for this spacetime in co-moving coordinates. 
The only non-zero component of the acceleration turns out to be 

ill = --. (23) 

X 

The non-vanishing rotation component becomes 

cuoi = 2-(-)^. (24) 

X Xq 
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The shear component is 

(7oi = -uoi (25) 
while the shear invariant becomes 

a = -4— = -iiil (26) 

With Raychaudhuri's equation, the rate of change of expansion takes the 
following form: 

de A(A + 25 - 1), X ,2R 24 

= 4 ^ ^(_)2B-2A_ 27 

ar x"^ x"^ Xq 

For the self-similar solution given by Eq.(9), the expansion scalar is given 

by 

e = i^. (28) 

The acceleration and the rotation turn out to be zero for this metric. The 
shear components are given as 

(2B-,A)U a^='-^-*''\ir = ,^, (29) 
and the shear invariant becomes 

+ 3A^ - AAB , , 
■ (30) 



^2 

The rate of expansion takes the form 
dB 1 



— ^(7^2 + 105^ + 3^ + 65-8^5). (31) 
or ?)t^ 



B. Class 2 



Class 2 has five solutions that are given by Eqs.(7), (10), (12), (14) and (16). 
In this class, for the first solution, the expansion scalar and rotation are zero. 
The only acceleration component is given by 

iix^ . (32) 

X 
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The only component of the shear becomes 



X Xo 

and the shear invariant takes the form 



aoi = -2-{-f, (33) 



a = -4-. (34) 



The expansion rate takes the following form: 

^ = ;.(^ + i). (35) 

For the second solution, given by Eq.(lO), the expansion scalar is given 

by 

e = ^. (36) 

The acceleration and the rotation are zero. The non-zero components of the 
shear are ^ 

C^ll = - — ( — )^'^ = (^22 = (733, (37) 
t to 



and the shear invariant becomes 



.= ^. (38) 
The rate of change of expansion takes the form 

- = -^(A + l). (39) 

The third solution in class 2 has the following expansion scalar: 

e = ^. (40) 

ax 

The only acceleration component is given by 

iii^- (41) 

X 
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while the rotation turns out to be zero. The non-vanishing components of 
the shear are 

(y\\ = — > cr22 = e " = 0-33, (42) 

iax iaxQ 



and the shear invariant becomes 



(43) 



The rate of change of expansion has the following form: 

For the fourth solution, given by Eq.(14), the expansion scalar is given 

by 

e = ^. (45) 

The acceleration and the rotation are zero for this metric. The components 
of the shear are 

(722 = « = (733 (46) 

while the shear invariant becomes 

The expansion rate is given, by using Raychaudhuri's equation, as 

53 

The last solution in this class, given by Eq.(16), has the expansion scalar 

e = 2a. (49) 

The acceleration and the rotation become zero. The components of the shear 
are 

(T22 = -ye2'^^*-^^)=a33, (50) 
and the shear invariant becomes 

o-^. (51) 

The expansion rate is 

de_ 26a^ 

d7 " ""9~- ^^^^ 
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2. Kinematic Properties in Non-co- moving Coordinates 

This section is devoted to discussing the same properties of the self-similar 
solutions in non-co-moving coordinates only in the radial direction. 

A. Class 1 

For the first metric in class 1, given by Eq.(6), we obtain non-zero expansion 
as follows: 

e = JA±lEl. (53) 

X 

The non-zero components of the acceleration turn out to be 

Uo=-{-f, iii^--. (54) 

X Xf) X 

The non-zero rotation component is 

oJoi = -(— )^ = iio. (55) 

X Xq 

The components of the shear are 

X Xq 3x Xq 

an = ^{A-B), a22 = -^(-)'^(A + 85) = a33, (56) 
6x 6x Xo 

and the measure of anisotropy is given by 

(7 = (4M2 + 16 AB + 70^2) . (57) 

The rate of change of expansion using Raychaudhuri's equation becomes 

^ = — ^(101^2 + 1345^ + 62AB - 18B) - -. (58) 
dr 9x^ X 

For the second metric in this class, we obtain non-zero expansion as fol- 
lows: 

e^^. (59) 
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The non-zero components of acceleration are 

A A,t 



^o^-j, = (60) 



and the non-zero rotation component is 



u;oi = ji^J^ = ^1. (61) 

The non- vanishing components of the shear are 

_ 2A _ {1QA + 2B) t ^ 

{10A-4B) t {A -AB) t 

— 3^ — y ' ^-=^^y =^^^' ^^'^ 

and the shear invariant is given by 

2 

a = — 
9t 

The expansion rate becomes 



2 (19^2 - 68AB + 225^). (63) 



^ = -H^(41^' + 385' - 106AB + 185) - -. (64) 

B. Class 2 

For the first metric in this class, we obtain non-zero expansion as follows: 

e = -— . (65) 

X 

The non-zero components of the acceleration turn out to be 

uo^-{-f: ^1 = --, (66) 

X Xq X 

and the non-zero rotation component is 

cooi = -(-)^ = ^0. (67) 

X Xq 
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The non- vanishing components of the shear are 

(Too = 4 — ( j , (7oi = ( j , (722 = ( ) = (T'as 

X Xq X Xq X Xq 

while the measure of anisotropy is given by 

30^2 

^ = — ~- 

x^ 

The rate of change of expansion becomes 

dQ A,^^^ A 
— = 33^ + 2 --. 
or x'^ X 

The second solution in this class gives the non-zero expansion 

t 

The non-zero components of the acceleration are 

A. . A. t .A 

and the non-zero rotation component is 

The components of the shear become 

Coo — — —■, ctqi — ~r\~r) ■> 

t T lq 

, t A t 

(711 = T^t^^ ' '^'^'^^ ~J^t^' =^33, 

and the shear invariant is given by 

The rate of expansion becomes 

de A A, ^ 

d7 = -T + ?P^-2)- 
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For the third solution, the non-zero expansion is given by 

X ax 

The components of the acceleration are 

1 1 

iio = — , -iii = — , 

Xq X 

and the rotation component is 

1 

Xo 

The components of the shear take the form 

Ax Axo 2 

c^oo = — cm = - — , ctqi = (3a + xq), 

Xo oax 6axo 

X . s 2-^ 



c^22 = ^(9a + 4xo)e » = (733 



Saxo 

while the measure of anisotropy is 

(7 = ^J^(135a2 + 22x0^ + eOaxo). 
The expansion rate becomes 

? = -7r^(261a' - - 156axo) - - + ^^(a' - Xo^)e' 

The fourth self-similar solution (Eq.(14)) gives the expansion as 

^ _ 3a - 2to 
at 

The acceleration components are 



t to 

while the non-zero rotation component is 

1 

UJoi = -— = -Ui. 

to 
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We obtain the following non-vanishing components of shear: 



2t 2 
c^oo = -T) CFii^ —-——-^{2,a + 2tQ), (7oi = — (a + io), 
t 6ato ato 



and the shear invariant is 

1 



The expansion rate becomes 

d& 1 2 



(^22 = -FV^^^" (3« + 8^o) = CTss, (86) 

OCLLq 



(45a2 + 44i2 + 24aio). (87) 



dr t daH"^ 



(45a^ + 19t^-6aio). (88) 



The expansion, the acceleration, and the rotation vanish for the last so- 
lution in this class. The shear components are 

(722 = -4ae'"(*+^) = (733, (89) 

and the shear invariant is 

a = 32al (90) 

The rate of expansion becomes 

^ = -30a\ (91) 
dr 

IV. Singularities 

In this section, wc shall explore the singularities of the self-similar solutions 
in classes 1 and 2. The Kretschmann scalar is defined by 

K = RaicdR"'"', (92) 

where Rabcd is the Riemann tensor. For the solution given by Eq.(6), Eq.(92) 
reduces to 

K = —AA^{A^ + 1 + 2^2 - 2A) + 5^(352 + 2 - AB)]. (93) 
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It is clear that K diverges ai x — 0. It follows that the solution is singular 

at a; = 0. 

For the solution given by Eq.(7), the Kretschmann scalar becomes 

2/12 

K^^{6A^ + 3-2A), (94) 

which shows that the spacetime singularity lies at x — 0. 
The solution given by Eq.(9) has the Kretschmann scalar 

K = \{A\A - if + 2B^{B - if + 2A^B'^ + B^). (95) 

V 

From here, it follows that the spacetime is singular at t = 0. 

For the metric given by Eq.(lO), the Kretschmann scalar reduces to 

K^^{2A^ + 1-2A). (96) 

V 

It is clear that K diverges at t = 0. Hence, the solution is singular at i = 0. 
For the solution given by Eq.(12), we obtain 

6 

K^-T-M'-^o'f. (97) 

which gives the spacetime singularity at a; = 0. The solution given by Eq.(14) 
has the following Kretschmann scalar: 

K^^^{tl-a^f. (98) 

It follows that the spacetime singularity lies at t = 0. 

For the last solution, given by Eq.(16), the Kretschmann scalar reduces 

to 

K = Aa^. (99) 
It turns out to be constant, which shows that this solution is singularity free. 



V. Conclusion 

Self-similar solutions in GR are very important, and discussing their physical 
features is interesting. Keeping this point in mind, we explored some kine- 
matic properties and the singularity features of such solutions representing 
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plane symmetric spacetime. We discussed the properties both in co-moving 
and non-co-moving coordinates (only in the radial direction). This provided 
a comparison of these properties in the two coordinates. We explored the 
acceleration, the expansion, the rotation, the shear, the rate of change of 
expansion, and finally the singularity. 

First, we discussed the class 1 which has two solutions. For the first solu- 
tion, given by Eq.(6), we found zero expansion in co-moving coordinates and 
positive/negative expansion in non-co-moving coordinates depending upon 
the values of the constants A and B. The acceleration and the rotation had 
only one component in co-moving coordinates while two components of the 
acceleration existed in non-co-moving coordinates. The rotation component 
remained the same in both coordinates, except for a factor of one-half in 
non-co-moving coordinates. We found the shear invariant to be negative in 
co-moving coordinates and positive/negative in non-co-moving coordinates, 
depending on the values of the constants. The rate of change of expan- 
sion could be positive/negative, depending upon the values of A and B in 
co-moving coordinates and negative in non-co-moving coordinates. For the 
second solution, given by Eq.(9), we obtained the same expansions in both 
coordinates, which could be positive/negative. These solutions had vanish- 
ing acceleration and rotation in co-moving coordinates and one rotation and 
two acceleration components in non-co-moving coordinates. The shear in- 
variant was positive in co-moving coordinates whereas the rate of expansion 
was negative. In non-co-moving coordinates, these quantities could be posi- 
tive/negative, depending on the values of constants. 

Now, we discuss the solutions of class 2. We notice that in co-moving 
coordinates, the solutions given by Eqs.(7) and (12) have non-zero acceler- 
ation component while the remaining solutions have zero acceleration. In 
non-co-moving coordinates, all solutions have non-zero acceleration compo- 
nents. Also, all solutions have zero rotation in co-moving coordinates while 
only one solution (Eq.(16)) has zero rotation in non-co-moving coordinates. 
The solution in Eq.(7) has zero expansion in co-moving coordinates whereas 
all other solutions have positive/ negative expansion. In non-co-moving coor- 
dinates, the expansion is positive/negative for the solutions in Eqs.(7), (10), 
(12), and (14) and zero for the solution in Eq.(16). 

The shear invariant is negative in co-moving coordinates and positive in 
non-co-moving coordinates for the solution in Eq.(7). The expansion rate 
is positive in co-moving coordinates, but it is negative in non-co-moving 
coordinates. For the solution in Eq.(lO), we have a positive shear scalar 
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in co-moving coordinates and a negative one in non-co-moving coordinates. 
The rate of expansion is negative in co-moving coordinates and positive in 
non-co-moving coordinates. For the solution in Eq.(12), the shear invariant 
is positive in both coordinates, and the expansion rate is negative in both 
coordinates. The shear scalar is positive and the rate of expansion is negative 
in both coordinates for the solution in Eq.(14). For the solution in Eq.(16), 
the shear invariant is positive, and the expansion rate is negative in both 
coordinates. 

Finally, we discuss the singularity structure for these solutions. The solu- 
tions given by Eqs.(6), (7), and (12), are singular at x — while the solutions 
given by Eqs.(9), (10), and (14) are singular ait — 0. The solution in Eq.(16) 
is singularity free. 

We have noticed from the above discussion that the kinematic quantities 
are relatively simple in co-moving coordinates as compared to those in non- 
co-moving coordinates. It is worth mentioning that the expansion of the 
solution in Eq.(16) turns out to be positive in co-moving coordinates, but it 
is zero in non-co-moving coordinates. 
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